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Introduction  and  Notation 


1-k 

Using  the  concept  of  aliasing,  Snider  [6]  obtains  an.  U(n  ) 
estimate  of  the  uniform  accuracy  of  the  n-point  trigonometric  inter- 
polants  of  periodic  functions  for  k > 2,  In^jroving  the 


•1/2  2 

0(n  ' ) estimate  for  C functions  presented  in  Isaacson  and  Keller 


[2].  Kreiss  and  Oliger  [k]  use  aliasing  to  show  that  if  the  Fourier 
coefficients  v(5)  of  a periodic  function  v(x)  satisfy 


•A 

v(5)  = U(  Ul~^)  with  p > 1,  then  the  trigonometric  interpolants  of 


V uniformly  approximate  v to  order  Cl(n  This  also  gives  an 

l*k 

cr(n  ) estimate  for  Cr  functions  since  the  largest  P we  can  use 
in  general  is  P = k.  We  use  aliasing  and  a different  property  of  the 
Fourier  coefficients  of  functions— the  fact  that  is  contained 

in  the  Sobolev  space  — to  obtain  an  o(n^/^”^)  estimate  for 
k > 1. 

2 

In  [5],  Kreiss  and  Oliger  estimate  the  L accuracy  of  trigonometric 
intei^olants  and  their  derivatives  for  functions  in  Sobolev  slices. 

This  ijaper  applies  their  approach  and  an  extension  of  a theorem  appearing  . 
in  Zygjnund  [7l  to  obtain  an  o(n^/^^”®)  estimate  of  the  uniform 


accuracy  of  the  m-th  derivatives  of  trigonometric  intei^olants  of  func- 


8 1 

tions  in  the  Sobolev  spaces  H for  s > ^ + m.  By  similar  methods 


we  obtain  an  o(n™  estimate  for  functions  in  whose  k-th 
derivatives  have  absolutely  converging  Fourier  series  if  k > m,  and 
we  show  that  these  two  estimates  are  sharp.  We  also  obtain  an 

estimate  for  functions  in  the  Holder  space  if 

0 < a < 1 and  k + a > ^ + m.  These  resvtlts  remain  valid  if  we  replace 
the  trigonometric  Interpolant  by  its  K-th  partial  sum,  replacing  n by 


4 


2 


L lv(|)|  < 

!=-« 

For  0 < a < 1,  let 

[v]  = sup 

x,y6]R 

For  an  Integer  k > 0,  is  the  set  of  functions  v(x)  e 

such  that  < ». 

If  V € A,  then  v is  equal  a.e.  to  a continuous  function.  Since 
we  are  interested  in  interpolation,  we  will  tacitly  assume  that 
A c C°  and  similarly  that  H®  c C°  for  s > ^.  For  an  integer 
k > 1,  ^ is  the  set  of  functions  v(x)  such  that  D^v  e and 
thus  c See  Agmon  [1]  for  a discussion  of  derivatives. 


2.  Trigonometric  Interpolation 

We  state  some  veil  known  results  on  trigonometric  interpolation. 
These  appear  in  this  form  for  odd  n in  Kreiss  and  Oliger  [4].  See 
also  Isaacson  and  Keller  [2]  and  Zygniund  [7]. 


A.  n is  odd.  Let  N > 0 be  an  Integer  and  h = Trrrr  and  let 
2N+1 

= vh  for  V = 0,1,2, ... ,2N.  There  is  a imique  trigonometric 
polynomial  Ifgv(x)  of  order  at  most  H which  interpolates  v(x)  at 


the  points  x^  for  0 < v < 2N  given  by 


2jti|x 


where 


^ |=-N 


2N  -2ni5x 

a(|)  = h r v(x  )e  '' 

v=0 


The  effect  called  aliasing  is  the  fact  that 


(3)  a(6)  - r v{5  + J(2N+1))  |s|  < N 

J=-0D 

provided  that  the  Fourier  series  for  v(x)  converges  at  the  points  x 

V 

for  0 < V < 2N. 

Following  the  notation  of  zygmund,  define  for  1 < K < N 


j,v(x)  = r a(S)( 
f=-K 


2«i5x 


where  a(5)  is  given  by  (2).  ^ IjV  is  the  K-th  partial  sum  of  I^v, 

and  Ijg  = Ijj'v*  If  v(x)  is  real-valued,  so  is  j^v. 
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B.  N Is  even.  Let  N > 0 be  an  integer  and  h = ^ and  let 

= vh  for  0 < V < 2N-1.  There  is  a unique  trigonometric  polynomial 
Ejjv(x)  of  order  at  most  N which  interpolates  v(x)  at  the  points 
x^  for  0 < V < 2N-1  given  by 


which  also  satisfies 


I 

S«v(x)  = r a(5)e' 
|=-N 


a(-N)  = a(N)  . 


2jti|x 


The  Z'  notation  indicates  that  the  first  and  last  terms  are  multiplied 
by  1/2.  The  coefficients  are  given  by 


2N-1 

a(0  = h Z v-(x  )e 
v=0 


Provided  that  the  Fourier  series  for  v(x)  converges  at  the  points 
X for  0 < V < 2N-1,  we  have 


a(g)  = E v(5  + J(2N)) 

J=-CB 


U1  <N 


Define  for  1 < K < N 


IV 

„v(x)  = E a(5)' 

5=.K 


2ni|x 


where  a(S)  is  given  by  (6),  and  let  Ejj 

real-valued,  so  is  Ej^  j^v  for  K < N.  If  w(x)  is  a trigonometric  poly- 
nomial of  order  at  most  N and  w(n)  = w(-N),  then  E^w  * w. 
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3.  Accuracy  Estimation 


» 


« 


Define 

6(v,m,N,K)  = |1d“v  - 
e(v,in,N,K)  = Hd^v  - d“( j^v ) ||^ 

The  m = 0 case  of  the  following  lemma  appears  in  Theorem  5.l6  of 
Chapter  10  in  Zygmund  [?]• 

Lemma  1.  Let  m > 0 be  an  integer,  and  suppose  that  u = D™v  e A. 
Then 

6(v,m,N,K)  < 2 E |u(|)l 

UI>K 


Proof.  Let 


(9) 


v^Cx) 


5=-K 


Vjj(x) 


_ -'z . X 2jti5x 

E v(  g )e 

UI>K 


(10) 


''l  ~ 


''r  = ^,K^R 


Then  V = and  = w^^  + Wj^.  Since  w^  = 

^ = ^r  - ''r 

so 

(12)  6(v,m,N,K)  < \\d\\\^  + I1d\|1^  . 

9/  (3), 

T 


1 


Vt,(x)  = E E V„(|  + J(2N+l))e^"^^'' 
« |=-K  J=-«  ^ 


K CO 

I|d\(L<  r l^nsT  E ly?  + J(2N+1))( 

|=-K  j=-co  ^ 


K CO 

< E E |2Jt(5  + J(2N+1))|“1vr(5  + J(2N+1))1 
5=-K  j=-»  ^ 


< E i2«iriyoi 


(13)  IId\IL<  e |u(i)| 

111  > K 


Also 


(1*^)  11d\1I„  < , e 1u(01 

III  > K 


Combining  (12),  (15),  and  (14)  gives  the  lemma. 


Lemma  2.  Let  m > 0 be  an  integer,  and  suppose  that  u = d'^v  e A. 
Then 


e(v,m,N,K)  < 2 E |u(|)|  for  K<N 

Ul  > K 

e(v,m,N,N)  < 2 E |u(0| 

111  > N 


Proof*  For  K < N,  the  proof  is  the  same  as  in  Lemma  1. 


i 


Using  (9)  with  K = N - 1 and  replacing  (10)  by 


“l  = Vl 


”e  ' Vr 


we  obtain 


(16) 


By  (7), 


e(v,in,N,N)  < ||d%||^  + ||d“wj^„^ 


Vj^(x) 


r'  S V (I  + j(2N))e^’'^^'' 
5=-N  J=-«o  ^ 


N 


< 2)  ' r )2n(5  + j(2N))riy  5 + J(2N))) 


5=-N  J=-® 


r |2«*J1V5)I 

I 


and  the  lemma  follows  as  in  the  proof  of  Lemma  1. 

Theorem  1.  Let  m > 0 be  an  integer  and  v e with  s > i 
Then  for  each  K, 


(17) 


where 


sup  6(v,m,N,K)  < CRj^( v )ic‘'/^ 


N > K 


2 (2n)“"® 
Vs  - ^ - m 


i 


+ m. 


9 


! 


, j 


and 


nJy)  = (,  E |2n5|^^lv(|)l^)V2  . 

UI>K 


Also 


(18) 


and 


(19) 


Note  that 


sup  e(v,m,N,K)  < CRj^(v)K^/^ 


-tm-s 


N > K 


e(v,m,K,K)  < CI^_^(v)  (K-1)^^^ 


since  v e E^,  ^ ss  K ->  ». 


Proof.  By  Lemma  1,  for  N > K we  have 


6(v,m,N,K)  < 2 L l2n«riv(0l 
|5|>K 


<2(  r |2«e|^^|v(0|^)^'^^(  L 

UI>K  U|>K 

^+2(m-s)  1/2 

< 2 Vv)(2n)-^(2 


and  (1?)  follows.  (l8)  and  (19)  follow  similarly  from  Lemma  2. 

Jr 

Theorem  2.  Let  k > m > 0 be  integers,  and  suppose  D v e A.  Then 
for  each  K, 
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riiMMMiiiiiiii 


(20) 


where 


and 


sup  5(v,in,N,K)  < Crj^(v)l^^“^ 


H > K 


m-k 


Also 


(21) 


and 


(22) 


C = 2(2n) 


" , 1^  |2«ir|v(|)l 

kl>K 


sup  €(v,m,N,K)  < Crj^(v)K™  ^ 


N > K 


e(v,m,K,K)  < Crj^_^(v)K®' 


Note  that  since  D v 6 A,  ^k(''^)  ^ K -♦  «. 


Proof.  By  Lemma  1,  for  N > K we  have 


6(v,m,N,K)<2  T |2«g|“lv(0l 

lll>K 


vffl-k 


< 2(2jtK )“*■'"  E l2ngl^lv(0l 

lll>K 


and  (20)  follows.  (2l)  and  (22)  follow  similarly  from  Lemma  2. 
Theorem  3.  Let  m >0  be  an  integer  and  v e with 


+ a > ^ + m.  Then  for  each  K, 


11 


J 


(23) 


sup  6(v,m,N,K)  < 
N > K 


where 


2a+i/2^-k 

^ gl/^4in-k-0£ 


Also 

(24)  sup  e(v,in,N,K)  < C[D^v]  ^ ^ 

N > K 


L 


k- 


Proof.  The  method  of  proof  is  similar  to  that  of  Bernstein's  theorem 
that  c A for  a > ^.  See  Katznelson  [3].  Let  u = D®v  and 

f = D^v.  If  t = j 2"'^  and  s''  < U 1 < 2'^'^^  then 
so  since 

f(x+t)  - f(x)  = L - l)f(5)e^"^^'‘ 

5=-co 

Parseval's  relation  implies  that 


2iti5t 


irif(5)r 


<5  llf(x+t)  - f(x)|l| 
< j llf(x+t)  - f(x)ll^ 


Ety  the  Schwarz  inequality. 


a''  < 1 6 1 < 2 


2''<  U|<2' 


2'^  < 1 5 I < 2'^’*'^  i2„g|2(k-m) 


< (2nf-^  2''(V2-H>i-k)(2 


2'’<  llf<2' 


< (2«f’^  2''^V2-^m-k-a)f^j 


Given  K,  let  j satisfy  2*^  < K < 2*^^^.  Then  hy  Lemma  1,  for 
N > K we  have 


6(v,m,N,K)<2  T |u(|)l 

UI>K 


“ v^j  2''<  uf<2'''^^ 


|u(Ol 


<2(2«r^[fl  I 2v(l/2-hn-k-a) 


Tn.v  ,pJxl/2-Hn-k-a 

< 2(2Jt)  ^ _ 2l/2-Hn-k-a 


and  (25)  follows  since  j > 2*^  and  ^ + m-  k-  a<0.  (2k)  follows 


similarly  from  Lemma  2. 
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4, 


. Sharpness  of  Estimates 

s 1 

Theorem  1 shows  that  if  v e H and  s > ^ + m,  then 

6(v,m,N,K)  and  e(v,m,N,K)  are  independent  of  N > K. 

k 

Theorem  2 shows  that  if  D v e A and  k > m,  then  6(v,m,N,K)  and 
e(v,m,N,K)  are  0(1^^“^),  independent  of  N > K.  We  prove  in  this 
section  that  these  estimates  are  shaip:  they  cannot  be  liqproved  for 
these  two  classes  of  functions. 


Theorem  4.  Let  {7^}  be  a sequence  of  positive  numbers  converging  to 
0.  Let  m > 0 be  an  integer,  and  s > ^ + m.  Then  there  exists  a 
v e such  that 


(25) 


lim  sup 
K — ♦ ® 


jj  > K 


CD 


Proof.  Let  p = 1 and  define  a strictly  increasing  sequence  {p.} 
of  positive  integers  inductively  such  that  for  J > 1#  if  J is  odd 
Pj  = 2pj  and  if  J is  even  Pj  is  a power  of  2 such  that 

(26)  7^  < 2“'^  for  V > Pj/^  • 

Define  the  sequence  {b^}  for  v > 1 hy 

4D  4D  OD 

Then  L h^  = E E b^  * 2 2"'^  < •. 
w=l  ''  Pj  $'’<Pj+^  J=0 
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Note  that  for  v > 1 since  Pj  > for  J > 0.  Let 

(28)  v(x)  = E (-1)''— i b 

v=l  (2nv)® 

OP  OP 

Since  E |2n|l^®lv(5)|^  = E b^  < <»,  v e H®.  Define  v_,  v_,  w^, 
5=-«>  v=l  L R L 

and  Wp  as  in  (9)  and  (10).  By  (U), 

(29)  «(v,m,N,K)  > ||d\1|^  - ||d\||^  . 

Now 

lD%(|)|=l.  £ (2«i5)“v(5)e’'^^l  = Z (2nv)“-«b 

UI>K  v>K  '* 

so 

(30)  llD%IL>  S (2nv)“"®b  . 

v>K  '* 

By  (5), 

Vj,(x)  = E a(Oe^’'^^'' 

^ |=-K 

where  for  | $ | < K, 

00  00 

a(5)  = E v„(5  + J(2N+1))  = E v(5  + j(2N+l)) 

J=—  “ J=1 

Since  2N  + 1 is  odd,  this  last  series  is  an  alternating  series  of 
terms  decreasing  in  absolute  value,  so 

15 


Hence 


la(0|  < lv(|  + 2N  + 1)1  . 


I|d\IL<  S |2«5ria(5)l 


< r |2rt(|  + 2N  + l)|®|v(|  + 2N  + 1)1 
|=-K 


2N+1+K 


E (2nv)”‘\ 


v=2N+l-K 


< r (2jrv)“"®b 


v=K+l 


since  the  form  a decreasing  sequence.  Combining  this  vith  (29) 

and  (50)  yields 


6(v,m,N,K)  > E (2nv)“‘®b  . 

~ v=3K+2  '' 


For  even  J > let  Kj  = Pj/^*  Then  since  Pj+.2  “ 2pj, 


6(v,m,N,K.)>  E (2rtv)®‘®b 


> S (2.>vr»(p,2J)-V2 

"j  Sv  <Pj„l 


p.  X 

J 


* 


- 


2 *-^-1 

1-p 

for 

P/^1 

log  2 

for 

P = 1 

so  If 


«(V,.,N,K  ) > 0^.^  2-V2(2„)”-V1/2«-« 


s-m  J 


Thus  (26)  iii5)lies  that 


-^75«rs  > 

Kj  J 


and  the  theorem  follows. 


Theorem  Let  {7^}  be  a sequence  of  positive  numbers  converging  to 
0.  Let  k > m > 0 be  integers.  Then  there  exists  a v with 


D V e A such  that 


(31) 


inf  6(v,m,N,K) 
n>K 

lim  sup  — ^ — 

K-.-  rj 


Proof.  Same  as  the  proof  of  Theorem  k with  the  following  alterations. 
Replace  s by  k throughout  the  proof.  Replace  (26)  by 


(26*) 


7 < 2 
V — 


-2J 


for  V > Pj/4 
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Define  b = 


for 


p.  < V < p 


J+1 


Then  E b < • and  E |2«5 |v(  5 ) | < • so  D^v  e A.  We  have 


v=l 


g=-00 


for  even  J 


«(v,m,N,K,)>  E (2nv)“’^ 
•>  - v-P^ 


> r 


(2nv)“"^Pj2'J)‘^ 


> (pj2*J)"-^(2nf'^  J 


.2P. 


r J 

c-m 


Pj  X 


= c^  2-j(2xr-v-^ 

k-m  ' ' 


— d.  2”*^K^”^ 
2 ^ 


Thus  (26')  lnplles  that 


6(v,m,N,K.) 

^ I 


■Kp 


I?-'' 


and  the  theorem  follows. 


The  following  lenna  Is  geometrically  obvious. 


Lenna^J.  Let  {3^)  be  a decreasing  sequence  of  positive  numbers 


converging  to  0,  Then  E 3 e 


,2nlv/5 


-T  V 

v=l 


converges  and 


‘ \l 


i 

! i 


1 i 
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Theorem  6.  Let  {7^}  be  a sequence  of  positive  numbers  converging 
to  0.  Let  m > 0 be  an  integer,  and  s > ^ + m.  Then  there  exists 
a V € such  that 


(32) 

and 

(33) 


inf  e(v,m,N,K) 

ii„  ,„p  

K 


ito  sup  ■ - . 


H . rj, 


If  k is  an  integer  vith  k > m,  then  there  exists  a v vlth  D^v  e A 
such  that 


(34) 


lim  sup 
K-*co 


inf  e(v,m,N,K) 

N>K.5i'N 


-k 


and 


(35) 


Urn  sup  . 


Proof.  The  proof  of  (32)  is  the  same  as  the  proof  of  Theorem  4 vlth 
the  following  alterations.  Replace  (28)  by 


v(x)  = t ^ b S®”*"' 

u=l  (2»v)“ 


For  N > K,  ve  have 


€(v,m,N,K)  > ||d\11.  - h\\\. 


L9 


where  is  given  by  (9)  and  Wj^  = Now 

|i>“v„(|)l  = 1 r (2«isf?(s)e'*””/’l  - r (2«v)”-\ 

UI>K  v>K  '* 


P\L>  r (2«v)“-\ 

v>K  ' 


By  (7), 


Wp(x)  = Z a(5)e 


2n±ix 


where  for  | 5 | < K, 


a(e)  = E vj5  + J(2N))  = Ev(5+J(2N))  . 

J=—  ^ J=1 


Suppose  5 I N.  Then  J(2N)  cycles  through  the  equivalence  classes 


mod  3f  so  by  Lemma  3, 


|a(|)|  < lv(S  + 2N)1  . 


Hence,  as  before. 


wi  a 


V=K+1 


and  the  rest  of  the  proof  goes  through,  establishing  (32). 

To  prove  (33)  for  this  v,  imitate  the  proof  of  Theorem  4 as  above 
with  the  following  changes.  Define  Vj^  and  v^^  by  (9)  with  K = N - 1, 
and  define  Wj^  and  w^  by  (15).  Then 

€(v,m,N,»)  > 1|D%1|,  - l|D%jil|,  . 
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f-,  ^>:  ifc«Nr»'-*r^.>^v  'i, 


As  above. 


\\^\L  > 2 (2«v)“’®b 


v>N 


^ (7), 


N 


where 


vJx)  = E'  a(|)e^"^^* 
^ |=-N 


a(£)  = L v(5  + J(2N)) 
j=l 

0» 

a(-N)  = a(N)  = E v(N  + J(2N)) 
J=0 


for  Ul  < N 


For  N * Kj  for  even  J > 4,  3 | N,  so  by  Lemma  3, 

|a(5)|  < |v(5  + 2N)|  for  |||  < N 

|a(-N)|  = |a(N)l  < |v(n)1  . 


Hence 


I1d\1L<  r l2niria(Ol 

5=-N 

N-1 

< E l2n(6  + 2N)riv(5  + 2N)1  + l2nNl®lv(N)l 
|=-N+1 

5N-1 

= E (2rrv)®-\ 
v=N  '' 


So 


€(v,m,N,N)  > E (2jtvr"®b 
” v=5N 


and  (33)  follows. 
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1 


(54)  and  (55)  follow  by  similar  alterations  to  the  proof  of 
Theorem  5» 


Remarks.  Theorem  4 shows  that  the  o(K^/^^”®)  estimate  of 
6(v,m,N,K)  given  by  Theorem  1 is  sharp  by  showing  that  there  is  no 
function  g(K)  going  to  0 faster  than  for  vhich 

6(v,m,N,K)  = Cfl(g(K))  for  all  v e . Note  that  we  can  obtain  a 
real -valued  function  in  satisfying  (25);  since  the  trigonometric 

interpolants  of  real-valued  functions  are  real-valued,  at  least  one 
of  the  real  or  imaginary  parts  of  the  v constructed  must  also  satisfy 
(25).  Similar  statements  hold  for  Theorem  5 and  6.  Also,  many  of  the 


details  of  the  constructions  are  for  convenience,  e.g.  making  the  Pj's 


powers  of  2,  and  placing  the  singularities  at  x = ^ in  the  odd  case  and 
2 


at  X = 


in  the  even  case. 
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resiilts.  For  exan^jle,  if  D v is  not  necessarily  continuous  but 


is  of  bounded  variation,  then  v(|)  = so 

•k  k*l 

l|v  - = CT(n  ).  Or,  if  D v is  absolutely  continuous  (or 

equivalently  if  D^v  e L^),  then  v( | ) = o(  and  Kreiss  and 

Oliger's  proof  shows  that  [jv  - = o(n^”^)  if  k > 1.  See 

Katznelson  [5]  and  Zygmund  [7]  for  discussions  of  the  growth  of 
Fourier  coefficients.  We  conclude  with  a table  of  estimates. 


If  D^V  € 

then  l|v  - w^ll^  = 

for 

o(n^"^) 

k > 2 

l" 

o(nV2-k) 

k > 1 

cO,0( 

Cr(nV2-k-«) 

k + a > ^ 

h" 

o(nV2-k-) 

k + s > j 

9 

B.V. 

k > 1 

• 

A 

o(n”^) 

k > 0 . 
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